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An AdS/CFT calculation of screening in a hot wind 
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One of the challenges in relating experimental measurements of the suppression in the number of 
J/xj) mesons produced in heavy ion collisions to lattice QCD calculations is that whereas the lattice 
calculations treat J /if) mesons at rest, in a heavy ion collision a cc pair can have a significant velocity 
with respect to the hot fluid produced in the collision. The putative J/ip finds itself in a hot wind. 
We present the first rigorous non-perturbative calculation of the consequences of a wind velocity v 
on the screening length L s for a heavy quark-antiquark pair in hot N = 4 supersymmetric QCD. 
We find L„(v,T) = f(v)[l -« 2 ] 1/4 /ttT with f(v) only mildly dependent on v and the wind direction. 
This L„(v,T) ~ L B (0,T)/ % /7 velocity scaling, if realized in QCD, provides a significant additional 
source of J/$ suppression at transverse momenta which are high but within experimental reach. 



Twenty years ago, Matsui and Satz suggested that 
because the attraction between a quark and an anti- 
quark is screened in a deconfined quark-gluon plasma, 
the production of J/ip mesons should be suppressed in 
sufficiently energetic nucleus-nucleus collisions relative to 
that in proton-proton or proton-nucleus collisions, since 
the screened interaction between a c and c immersed in a 
quark-gluon plasma would not bind them [f[. fn the in- 
tervening years, marked progress on many fronts has not 
changed this basic qualitative picture. On the experi- 
mental side, we now have data from the NA50 and NA60 
experiments at the CERN SPS and from the PHENIX 
experiment at RHIC that demonstrate the existence of a 
suppression 0. On the theoretical side, we now have ab 
initio calculations of the temperature-dependent poten- 
tial between a color singlet heavy quark and antiquark 
separated by a distance L |j| . This potential is as at 
T = for small L, but begins to weaken for L larger 
than some L s and flattens at larger L. The potentials 
obtained in these lattice calculations can be crudely char- 
acterized as indicating ~ 0.5/T in hot QCD with two 
flavors of light quarks 0] and L s ~ 0.7/T in hot QCD 
with no dynamical quarks Q. Furthermore, lattice QCD 
calculations of the Minkowski space J/ip spectral func- 
tion itself have now been done in quenched QCD [|| , and 
early results in QCD with dynamical quarks have also 
been reported These studies indicate that the J/ip 
meson ceases to exist as a bound state above a tempera- 
ture somewhere between 1.5 T c and 2.5 T c , in agreement 
with conclusions drawn based upon the screening poten- 
tial between static quarks Q ■ 

The multifaceted challenge, now, is to make quanti- 
tative contact between the lattice calculations and data 
from heavy ion collisions. One significant difficulty is 
that the lattice calculations treat a quark-antiquark pair 
in the quark-gluon plasma rest frame, whereas in a heavy 
ion collision a cc pair is not produced at rest. This chal- 
lenge becomes more acute in higher energy collisions: 
at LHC energies, cc pairs which if produced in vacuum 
would yield J/ip mesons with transverse momenta many 



times their rest mass will be copious. Even in collisions 
at SPS and RHIC energies, the collective flow developed 
by the hot medium in which the cc pair finds itself is con- 
siderable. A rigorous determination of the ^-dependence 
of the screening length L S (T) for a heavy quark anti- 
quark pair in a quark-gluon plasma moving with velocity 
v would therefore be a significant advance. We provide 
one, albeit for hot N = 4 super Yang-Mills theory. 

Af = 4 super Yang-Mills (SYM) theory is a confor- 
mally invariant theory with two parameters: the rank 
of the gauge group N c and the 't Hooft coupling A = 
5y M A c . We shall define the screening length L s be- 
low, based upon an analysis of a fundamental Wilson 
loop describing the dynamics of a color-singlet quark- 
antiquark "dipolc" moving with velocity v along, say, the 
^-direction through the hot strongly interacting JV = 4 
SYM plasma. In the rest frame of the dipole, which sees 
a hot wind blowing in the cc3-direction, the contour C of 
the Wilson loop is given by a rectangle with large ex- 
tension T in the ^-direction, and short sides of length L 
along some spatial direction. Evaluating this Wilson loop 
(whether ultimately in QCD or in Af — 4 SYM as we do 
here) will teach us about the L-dcpendent color singlet 
quark-antiquark potential and hence allow us to define a 
screening length in the presence of a hot wind. 

According to the AdS/CFT correspondence [1], in 
the large- iV c and large-A limits the thermal expectation 
value (W F {C)) for the Wilson loop in the absence of a 
wind velocity can be calculated using the metric for a 
5-dimcnsional curved space-time describing a black hole 
in anti-deSitter (AdS) space Calling the fifth dimen- 
sion r, the black hole horizon is at some r = r$ and 
we add a probe D3-brane extended along the x\, X2, X3 
directions at some r = A ro- The external quarks de- 
scribed by the Wilson loop are open strings ending on the 
probe brane. The prescription for evaluating (W F (C)} is 
that we must find the extremal action surface in the five- 
dimensional AdS spacetime whose boundary at r = A 
is the contour C in Minkowski space R 3 ' 1 . (W F (C)) is 
then given by expos' (C)] = exp[iE(C)T], with S the ac- 
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tion of the extremal surface Q. For the time-like Wil- 
son loop we analyze, S(C) is proportional to the time 
T, meaning that E(C) can be interpreted as the energy 
of the dipole. In the limit of infinitely heavy quarks, i.e. 
A — * oo, E(C) oc A but this divergence comes from the L- 
independent self-energy of the quark or antiquark taken 
separately. We are only interested in the L-dependent 
part of E(C), which is finite in the A — > oo limit, and we 
take this limit henceforth [l6| ■ 

To describe a hot wind in the a^-direction, we boost 
the five-dimensional AdS black hole metric, obtaining 



ds 2 = -Adt 2 + 2Bdtdx 3 



-Cdxl- 



l (d Xl 



-dxD + f^dr 2 
(1) 

where / = 4w ( 1 — ^ ) with R the curvature radius of 



(2) 



the AdS space and where we have defined 



r 



r 2 R 2> 



B = 



cosh 



n 



d 1 

*o7 



' 1'2 

and r| 



C = 



r 

R 2 
I sinh 



r 2 R 2 



4 2 2 
V = r 0I V > 



and rf 

with v the velocity of the wind, rj its rapidity, and 7 = 
l/Vl — v 2 . Here, the temperature T of the Yang-Mills 
theory at v = is given by the Hawking temperature of 
the black hole, T = -%r, and R and the string tension 

l/2ira' are related to the t'Hooft coupling by -^7- = \f\. 

The short side of C can be chosen to lie in the (xi,xs) 
plane, at an angle 9 € [0, 7r/2] relative to the ^-direction. 
is the angle between the dipole and the wind direction 
in the dipole rest frame. We parameterize the surface 
whose action S(C) is to be extremized by x M = x m (t, a), 
where a a = (t, a) denote world sheet coordinates. The 
Nambu-Goto action for the string world sheet is given by 
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with g a p = G^daX^d^x" the induced metric on the 
world sheet. This action is invariant under coordinate 
changes of er", and we can set r = t and a = x\ [l7j ] . 
Note that a takes values — ^sin0 < a < -jsinfl. Since 
T L, we can assume that the surface is translationally 
invariant along the r direction, i.e. ^(c, r) = x^(a). 
The Wilson loop lies at constant x-i, so X2(cr) = const. 
The boundary condition for X3 (a) can be taken as 
X3(±-| sin#) = ±-|cos#. For the bulk coordinate r(a), 
we implement the requirement that the world sheet has 
C as its boundary by imposing r (±4 sin 0) = A — > 00. It 
proves convenient to define y = r/r and z = x^Tq/R 2 = 
x^ttT and hence I = Ltq/R 2 = LttT and to define a 
rescalcd a — ar^/R 2 and promptly drop the tilde. The 
boundary conditions become 
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The action ([3]) now takes the form 
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FIG. 1: String world sheet for wind with velocity v = 0.7 
blowing at an angle 6 = 45° relative to the dipole. The solu- 
tion has integration constants p = 1.325 and q = 1.109, which 
correspond to 6 = 45° and £ = 0.689. (This I is the maximum 
possible for this v and 9.) a oc x\ extends from — (^/2)sin# 
to +(£/2)sin#. (a) y(cr). (b) z(a) — a is the deviation of 
the string world sheet away from z = a, the straight line at 
9 — 45° between the quark and the antiquark. 

with the constant prefactor given by K = yXT and with 
the Lagrangian 



C = 



cosh 77) ( 1 + 



y 



1 



+ z' 2 (y 4 -l), (6) 



where henceforth 



Jjr. The equations of motion ob- 
tained from this Lagrangian can be cast in the form 
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(y i - cosh 2 V )(y A - 1 - P 2 ) ~ ?V - 1) , (7) 
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where p and q are integration constants. Since y' depends 
only on y and the boundary condition for y is symmetric 
under a <-> — a, we find y(cr) to be an even function 
of a: it descends (y' < 0) for — 1/2 < a < 0, turns 
around at a = where y' = 0, and then ascends. Since 
z' depends only on y and the boundary condition for z 
is antisymmetric in a — > — a, we find z(a) to be an odd 
function of a. In particular, z = at a = 0. 

Eqs. ([7]) and (jHJ) with the boundary conditions ^ 
can be integrated numerically and in Fig. 1 we present 
an example. Somewhat counterintuitivcly, the projection 
of the string world sheet onto the (xi, X3)-plane is not a 
straight line connecting the quark and the antiquark, but 
rather has a sinusoidal form. This behavior arises for all 
values of 9 except 8 = or n/2 (wind parallel or perpen- 
dicular to the dipole) for which the projection is indeed 
a straight line. Furthermore, we see that even though 
there is a wind blowing in the z-direction, y(a) is even 
and z(a) is odd for any angle 9, meaning that the string 
world sheet is not dragged at all by this wind. This con- 
clusion is antithetical to that for the world sheet of an 
isolated string ending on a single quark or antiquark, an- 
alyzed in Ref. 12| , and in qualitative agreement with the 
conclusion that "mesons feel no drag" reached in a differ- 
ent context in Ref. The 9 — n/2 case is particularly 



3 




FIG. 2: The dependence of the screening length on the wind 
velocity v. (a) I = LirT, given in Eq. |9j as a function of 
the integration constant q for six different values of v: 0, 
0.5, 0.7, 0.8, 0.9, 0.95 (top to bottom). We see the peak 
of this curve, £ m ax, dropping with v. All curves are for a 
wind blowing in the direction perpendicular to the dipole. 
(b) f(v) = = LaivT^y versus v for a wind blowing 

perpendicular to the dipole (lower curve) or parallel to the 
dipole (upper curve). 



simple to analyze because the integration constant p van- 
ishes 18[. With p = 0, z(cr) is constant, starkly making 
the point that the world sheet is not dragged by the wind 
blowing in the ^-direction. 

We shall find that the screening length behaves simi- 
larly for all values of 9. For purposes of illustration, we 
present the analysis for 9 = ir/2. With z(a) = constant, 
we need only analyze the shape of the world sheet in the 
r-direction, governed by ([7]) which we can write as y' = 
\\J(.V A ~ !)(y 4 ~ Vt) u P° n defining y\ = cosh 2 ?] + q 2 . 
The world sheet stretches from y = oo down to y = y c , 
with y'(Q) = and y(fi) = y c . The integration constant 

q can now be determined from the equation | 
which becomes 



So d(J ' 



1 



y c 4 )(y 4 -i) 



The energy can be written as 



y 4 — cosh rj 



E = K I dy 



(9) 



(10) 



and can be made finite by subtracting the self-energy of 
an isolated quark and antiquark. 

In Fig. 2a we plot £ as a function of q from Eq. for 
several values of cosh 2 r\. Analysis of Eq. §§§ shows that, 



for any cosh 2 rj, i decreases at large q like I m 1.198/ y/q. 
This corresponds to an energy E oc —K^/q oc —K/l oc 
— y/X/L, and thus describes a quark and antiquark sep- 
arated by a small distance L interacting via an attrac- 
tive Coulomb potential. Small values of L can also be 
achieved by choosing q — > 0, as I oc q at small q. These 
world sheets have E oc const + L 2 , meaning that these 
small L solutions have much higher energy than those 
describing Coulomb attraction, and so are not of inter- 
est. Since I is everywhere positive and goes to zero for 
both small and large q, at some q m it has a maximum 
value £ max as illustrated in Fig. 2a. For I > £ max , no ex- 
tremal world sheet bounded by the Wilson loop C exists, 
with the exception of disconnected world sheets "hanging 
from" the quark or the antiquark Wilson line separately, 
describing their self energies. Hence, for I > ^ max there 
is no L-dependent potential between the quark and an- 
tiquark [19(. We can therefore define a screening length 
L s by L s = e max /irT Analysis of © shows that 

for v — 0, q m = 0.96 and £ max — 0.869 arc numbers 
of order 1 while for large cosh 2 77, q m w 1.42 cosh 77 and 
Zmax ~ 0.743/Vcosh»7 = 0.743(l-w 2 ) 1 / 4 . This motivates 
writing the screening length as 

i 8 = fd-^, 



(11) 



in so doing defining the function f(v). We plot f(v) in 
Fig. 2b, and find that its velocity dependence is mild 
meaning that the dominant w-dependence of L s is the 
factor (1 — w 2 ) 1 / 4 . For general 9, the results are very 
similar: the interaction is screened for L > L s where L s 
is given by pip with a slightly different f(v), shown in 
Fig. 2b for a wind blowing parallel to the dipole. As the 
angle 9 between wind and dipole changes from 0° to 90°, 
f(v) interpolates between the two curves in Fig. 2b. 

Our central result is that, in M = 4 SYM theory, 
the dominant dependence of the screening length of a 
dipole in a hot wind on the wind velocity is L s (v,T) ~ 
L s (0,T)/ v ^', with the remaining weak dependence de- 
scribed by the function f(v) in Fig. 2b. The dominant 
velocity dependence suggests that L s should be thought 
of as oc (energy density) -1 / 4 , since the energy density in- 
creases like 7 2 as the wind velocity is boosted. It turns 
out that L s (0,T) is within a factor of two of that for 
QCD. Although it would be interesting to sec whether 
/(0) is closer to that in QCD in more QCD-like theories 
with gravity duals, given the availability of reliable lat- 
tice calculations of the screening potential in QCD itself 
this is not a pressing issue. It would certainly be inter- 
esting to see whether L s (v,T) ~ L s (0,T)/ y/^y persists in 
other theories with gravity duals, as this would support 
its applicability to QCD [211 ]. 

If the velocity-scaling of L s that we have discovered 
holds for QCD, it will have qualitative consequences for 
quarkonium suppression in heavy ion collisions. For il- 
lustrative purposes, consider the explanation of the J/'fy 
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FIG. 3: A 1/^/7- velocity scaling of the screening length in 
QCD would imply a J/\t dissociation temperature Tdiss(p:r) 
that decreases significantly with pt , while that for the heavier 
T is affected less at a given pr- The T curve is schematic: 
we have increased Tdi ss (0) over that for the J/ip by a factor 
corresponding to its smaller size in vacuum. 

suppression seen at SPS and RHIC energies proposed in 
Refs. 0,111: lattice calculations of the gg-potential in- 
dicate that the J/^(1S) state dissociates at a temper- 
ature ~ 2.1T C whereas the excited Xc(2P) and t'(2S) 
states cannot survive above ~ 1.2T C ; so, if collisions 
at both the SPS and RHIC reach temperatures above 
1.2T C but not above 2.1T C , the experimental facts (com- 
parable anomalous suppression of J production at the 
SPS and RHIC) can be understood as the complete loss 
of the "secondary" J/^'s that would have arisen from 
the decays of the excited states, with no suppression at 
all of J/vff's that originate as J/*'s. Taking eq. (fTTj) 
at face value, the temperature Tdi SS needed to dissoci- 
ate the J/^ decreases oc (1 — w 2 ) 1 / 4 . As can be seen 
from Fig. [3J this indicates that J/\l/ suppression at RHIC 
will increase markedly (as the J/'I^IS) mesons them- 
selves dissociate) for J/^P's with transverse momentum 
Pt above some threshold that is at most ~ 9 GeV and 
would be ~ 5 GeV if the temperatures reached at RHIC 
are ~ 1.5T C [23]. These illustrative considerations point 
to a novel quarkonium suppression pattern at transverse 
momenta above 5 GeV, a regime that is within experi- 
mental reach of future high- luminosity runs at RHIC and 
that will be studied thoroughly at the LHC. If the tem- 
peratures reached at the LHC are, say, ~ 3T C , the LHC 
could discover T suppression, but only at high enough 
Pt- 
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